Create representation of a sequence. Begining and end points are determined for stretches of nucleotides. Stretches are determined by querying each nucleotides in a sequence to determine if it is represented in the database of characters (chars).
proc.chromR(x, win.size = 1000, verbose = TRUE)regex.win(x, max.win = 1000, regex = "[acgtwsmkrybdhv]")
seq2rects(x, chars = "acgtwsmkrybdhv", lower = TRUE)
var.win(x, win.size = 1000)
gt2popsum(x)
gt.to.popsum(x)
The function regex.win() is used to generate coordinates to define rectangles to represent regions of the chromosome containing called nucleotides (acgtwsmkrybdhv). It is then called a second time to generate coordinates to define rectangles to represent regions called as uncalled nucleotides (n, but not gaps).
The function gt2popsum is called to create summaries of the variant data.
The function var.win is called to create windowized summaries of the chromR object.